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Abstract 

Abundances for C^®0, CS, NH3, H2CO, C3H2, and 
N2H+ and upper limits for the abundances of and 

HCO’*' are derived for gas within the Bok globule B68 
using our own BIMA array data, single-dish data from 
the literature, and Monte Carlo radiative transfer mod¬ 
els. B68 has had its density structure well determined, 
removing a major uncertainty from abundance determi¬ 
nations. All abundances for B68 are lower than those de¬ 
rived for translucent and cold dense clouds, but perhaps 
only significantly for N2H+, NH3, and C3H2. CS deple¬ 
tion toward the extinction peak of B68 is hinted at by 
the large offset between the extinction peak and the po¬ 
sition of maximum CS line brightness. C^®0 and N2H+ 
abundances derived here are consistent with recently de¬ 
termined values at positions observed in common. 
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1. Introduction 

We examine the abundances of several molecular species 
in the Bok globule Barnard 68 (B68; Barnard 1919), lo¬ 
cated at 125 pc in Ophiuchus (de Geus, de Zeeuw, & Lub 
1989) with no evidence for internal star formation. Figure 
1 shows an i?-band image of B68 obtained from the Al- 
adin sky atlas, revealing its compact morphology. Using 
the extinction of the background Galactic bulge popu¬ 
lation, Alves, Lada, & Lada (2001; ALLOl) found a ra¬ 
dial column density profile for B68 that was well matched 
with that of an isothermal Bonnor-Ebert sphere, specifi¬ 
cally one with a near-critical center-to-edge density con¬ 
trast of 16.5 (Bonnor 1956; Ebert 1955). This robust den¬ 
sity determination removes a major uncertainty towards 
determining the molecular abundances in B68. Recently, 
Bergin et al. (2002) used single-dish maps of B68 in con¬ 
junction with models incorporating the ALLOl Bonnor- 


Ebert density structure to ascertain radial variations in 
C^®0 and N2H’*' abundance, with the lowest values oc¬ 
curring toward the highest extinctions. We explore fur¬ 
ther abundances of these and other species in B68, also 
incorporating the ALLOl density structure into models. 
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Figure 1. ESO MAMA R-band image of B68 from Aladin. 
The extinction peak is denoted by a diamond. The central 
pointing position of the BIMA data is shown as a cross 
and the primary beam FWHMs at 89.2 GHz and 110.2 
GHz are shown respectively as outer and inner dashed cir¬ 
cles. The positions and beam FWHMs of single-dish line 
observations from the literature are shown respectively as 
triangles and dotted circles. The “X” denotes the position 
of B68 from Glemens & Barvainis (1988). 
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2. Observations 

Table 1 summarizes the data used to determine molecu¬ 
lar abundances in B 68 , listing the lines used, the source 
observatory, the angular resolution, and the angular off¬ 
set of the observation from the extinction peak. Figure 1 
shows the position observed using the BIMA millimeter 
interferometer, approximately 4" north of the extinction 
peak. Figure 1 also shows the minimum and maximum 
FWHM sizes of the BIMA primary beams at the frequen¬ 
cies of 1-0 and HCO’*' 1-0. Line emission was not 

detected with the BIMA array at any observed frequency 
or position toward B 68 . Table 1 also lists studies where 
detections of lines were made toward B 68 from single¬ 
dish observations having angular resolutions <60^b Fig¬ 
ure 1 shows the resolutions of these data and the positions 
where line characteristics were reported. Note that these 
positions, typically the locations of peak line intensity, 
vary widely across B 68 . 

3. Abundances 

Abundances were estimated using the ID Monte Carlo 
code of Hogerheijde & van der Tak (2000) to calculate 
radiative transfer, molecular excitation, and line emission 
through a model of B 68 , and varying abundances until 
the output matched the observed data. Our B 68 model 
had the thermal and density structures from ALLOl, i.e., 
an isothermal sphere of T = 16 K with an outer radius 
of 12 500 AU and a “Bonnor-Ebert parameter” ^max of 
6.9 (i.e., center-to-edge density ratio of 16.5). The central 
density is ~2.5 x 10® cm“®. Constant radial abundances 
were assumed. Although abundance gradients within B 68 
are quite possible (see Bergin et ah), they should not 
produce large differences in the derived abundances, given 
the small sizes of the beam widths relative to the spatial 
extent of B 68 . The velocity field of the gas was assumed 
to have a turbulent line width of 0.4 km s“® (see Wang 
et al.) with no bulk motions. 

To compare models with interferometer data, “visi¬ 
bility” datasets were produced from the model cubes us¬ 
ing the baselines from the original data. These datasets 
were inverted, cleaned and restored in the same manner 
as the actual data, and velocity-integrated line intensities 
calculated. By properly accounting for spatial filtering, 
we estimated BIMA would have recovered only ~20% of 
the flux emitted from an ALLOl Bonnor-Ebert sphere in 
B 68 . To compare models with single-dish data, the model 
cubes were convolved with Gaussians of widths appropri¬ 
ately representing the resolutions of the respective obser¬ 
vations. 


Table 1. Summary of Line Observations of B68 


Line 

Site 

Beam 

(" X ") 

Offset 

(") 

Ref. 

1-0 

BIMA 

17.1 X 4.7 

4.03 

1 

1-0 

BIMA 

17.7 X 4.5 

4.03 

1 

2-1 

CSO 

30 X 30 

59.3 

2 

CS 2-1 

FCRAO 

46 X 46 

66.2 

3 

HCO+ 1-0 

BIMA 

18.5 X 7.0 

4.03 

1 

N 2 H+ 1-0 

BIMA 

17.7 X 6.6 

4.03 

1 

N 2 H+ 1-0 

Haystack 

18 X 18 

16.3 

4 

NH 3 (1,1) 

Effelsberg 

40 X 40 

54.5 

5 

H 2 CO 3i2-2ii 

CSO 

30 X 30 

59.3 

2 

C 3 H 2 2i2-1oi 

Haystack 

18 X 18 

16.3 

4 


References — (1) this work, (2) Wang et al. 1995, (3) Laun- 
hardt et al. 1998, (4) Benson, Casefli & Myers 1998, (5) 
Lemme et al. 1996. 


Table 2 lists the fractional abundances of the sam¬ 
pled molecular species at various positions toward B 68 . 
The BIMA data yield the first upper limits to the abun¬ 
dances of HCO"'' toward the extinction peak of B 68 , as 
well as upper limits to the abundances of ^®CO, C ^®0 
and N 2 H+ at the same location. In addition, single-dish 
data of B 68 from the literature yield abundance values 
of H 2 CO, CS, NH 3 , and C 3 H 2 toward B 68 for the first 
time, but at positions offset from the extinction peak. Fi¬ 
nally, other single-dish data of B 68 from the literature 
yield new abundance values of C ^®0 and N 2 H+ at some 
of these latter locations. 

Table 2 lists fractional abundances from the literature 
for the species studied here, in clouds similar in charac¬ 
ter to B 68 . The first set, for “translucent clouds,” are 
those listed by Turner (2000) for round Clemens & Bar- 
vainis clouds with edge-to-center visual extinctions of 2 . 0 . 
The second set, denoted for “cold dense clouds,” are those 
compiled by Ohishi, Irvine, & Kaifu (1992) toward TMC- 
1 or L134N. Turner notes that these abundances should 
not be regarded as universal as they can vary by an order 
of magnitude both between different translucent clouds 
and within larger, dense clouds. Every value of molecular 
abundance we derive for B 68 is less than the lowest value 
derived for other clouds. However, only the N 2 H+, NH 3 , 
and C 3 H 2 abundances differ by over an order of magni¬ 
tude, and so only these species may be arguably depleted 
in B 68 , at least at the locations observed. However, Fig¬ 
ure 1 and Table 2 show that less-discrepant abundances 
are found at positions relatively far from the extinction 
peak but more-discrepant abundances are found at posi- 
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Table 2. Abundances of Various Molecular Clouds 


Species 

B 68 

Translucent 

Clouds 

TMC-1 

L134N 

i*CO 

<1.3(-7) 


8.9(-7) 

8.9(-7) 

c®«o 

<1.3(-7) 


1.6(-7) 

1.6(-7) 

c®®o 

3.0(-8) 


1.6(-7) 

1.6(-7) 

CS 

4.0(-10) 

l.l(-9) 

l(- 8 ) 

l(-9) 

HCO+ 

<1.4(-10) 

2(-9) 

8(-9) 

8(-9) 

N 2 H+ 

<1.3(-10) 

l(-9) 

5(-10) 

5(-10) 

N 2 H+ 

2 . 0 (-ll) 

l(-9) 

5(-10) 

5(-10) 

NH 3 

7.0(-10) 

2.1 (- 8 ) 

2 (- 8 ) 

2(-7) 

H 2 CO 

4.0(-10) 

6.3(-9) 

2 (- 8 ) 

2 (- 8 ) 

C 3 H 2 

1 . 2 (- 11 ) 

3.6(-8) 

3(-8) 

2(-9) 


Note — Derived from abundance assuming = 

90 or 0/1*0 = 500. 

Note — a{—b) denotes a x 10“® 


tions closer to the extinction peak (except notably NH 3 ). 
This pattern suggests C 3 H 2 may be also depleted by some 
process related to extinction, e.g., the sublimation of gas- 
phase molecules onto grains, as suggested by Bergin et 
al. for Ci®0 and N 2 H+. This same idea may also explain 
how CS and H 2 CO appear relatively undepleted in the 
outer, less-extincted radii of B 68 . Moreover, the lack of 
bright emission in CS near the extinction peak hints that 
CS may be also depleted at high extinction in B 68 . 

For their abundance estimates for B 68 , Bergin et al. 
assumed a Zucconi et al. temperature profile. They found 
the Ci®0 abundance in B 68 rises from very low values at 
Ay « 1 to a peak of 1 x 10 “'" at Ay = 2 , and decreases 
to 1 X 10“® at Ay > 20, a contrast of ~100. Also, they 
found the N 2 H’'' abundance rises from very low values at 
Ay « 1 to a peak of 6 x 10“ii at Ay = 3, and decreases 
to 3 X 10“ii at Ay > 20, a contrast of ~2. Any evidence 
of N 2 H+ depletion is remarkable, given its oft-described 
utility as a non-depleting probe of dense core interiors 
(e.g., see Tafalla et al. 2002.) 

Adopting the same temperature profile in our models, 
derived abundances are increased by factors of 2-3 from 
those listed in Table 2. Our interferometer data of C^®0 
and N 2 H''' do not provide much additional support for the 
abundance models of Bergin et ah, except that the non¬ 
detection of compact line emission with BIMA suggests 
the low abundances found by Bergin et al. are not due 
to the beam dilution of small clumps of relatively abun¬ 
dant material. Our derived upper limits for C^®0 and 
N 2 H+ from the BIMA data are not particularly low, but 


they remain consistent with the still-lower values found 
by Bergin et al. throughout the core. 

The single-dish C^®0 and N 2 H+ data from the litera¬ 
ture only pertain to one line-of-sight per transition toward 
B 68 , but provide strong support for the abundances de¬ 
rived by Bergin et al., again assuming the same Zucconi 
et al. temperature profile. The C^®0 literature data, from 
a position 59/4 offset from the extinction peak, yield an 
apparent abundance of 6-9 x 10“®, quite consistent with 
their 3-9 x 10“® abundances at this position. Also, the 
N 2 H+ literature data, from a position 16/3 offset from the 
extinction peak, yield an apparent value of 4-6 x 10“^®, 
slightly larger than but still consistent with their 3-4 x 
10 “^^ abundance at this position. 

Acknowledgements 

JD was supported by the Radio Astronomy Laboratory. MRH 
was supported by the Miller Institute for Basic Research in 
Science. We thank Charles J. Lada, Joao Alves, Tracy Huard, 
and Jon Swift for valued discussions. This research made use 
of the SIMBAD database and Aladin, both operated by CDS, 
Strasbourg, France. 

References 

Alves, J. F., Lada, C. J., & Lada, E. A. 2001, Nature, 409, 
159 (ALLOl) 

Barnard, E. E. 1919, ApJ, 49, 1 

Benson, P. J., Caselli, P, & Myers, P. C. 1998, ApJ, 506, 743 
Bergin, E. A., Alves, J., Huard, T., & Lada, C. J. 2002, ApJ, 
570, LlOl 

Bonnor, W. B. 1956, MNRAS, 116, 351 
Clemens, D. P., & Barvainis, R. 1988, ApJS, 6 
de Geus, E. J., de Zeeuw, P. T., & Lub, J. 1989, A&A, 216, 
44 

Ebert, R. 1955, Z. Astrophys., 37, 217 

Hogerheijde, M. R. & van der Tak, F. F. S. 2000, A&A, 362, 
697 

Launhardt, R., Evans, N. J., H, Wang, Y., Clemens, D. P., 
Henning, Th., & Yun, J. L. 1998, ApJS, 119, 59 
Lemme, C., Wilson, T. L., Tieftrunk, A. R., & Henkel, C. 
1996, A&A, 312, 585 

Ohishi, M., Irvine, W. M., & Kaifu, N. 1992, in lAU Sympo¬ 
sium 150, Astrochemistry of Cosmic Phenomena: Proceed¬ 
ings of the 150th Symposium of the International Astro¬ 
nomical Union, ed. P. D. Singh (Dordrecht: Reidel), 171 
Tafalla, M., Myers, P. C., Caselli, P., Walmsley, C. M., & 
Comito, C. 2002, ApJ, 569, 815 
Turner, B. E. 2000, ApJ, 542, 837 

Wang, Y., Evans, N. J., H, Zhou, S., & Clemens, D. P. 1995, 
ApJ, 454, 217 

Zucconi, A., Walmsley, C. M., & Galli, D. 2001, A&A, 376, 
650 



